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Abstract. One of the fundamental astrobiology questions is how life has formed in our Solar 
System. In this context the formation and stability of abiotic organic molecules such as CH4, formic 
acid and amino acids, is important for understanding how organic material has formed and survived 
shocks and energetic particle impact from winds in the early Solar System. Shock waves have been 
suggested as a plausible scenario to create chondrules, small meteoritic components that have been 
completely molten by energetic events such as shocks and high velocity particle impacts. We study 
here the formation and destruction of certain gas-phase molecules such as methane and water during 
such shock events and compare the chemical timescales with the timescales for shocks arising from 
gravitational instabilities in a protosolar nebula. 
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INTRODUCTION 

Recent observations have shown that protoplanetary disks are ubiquitous in locations 
where new stars are forming. At the same time, astronomers found that at least 10% of 
nearby Solar like stars harbor planetary systems very different from our own, namely 
with Jupiter-like planets on Earth-like orbits. These findings have revived discussions 
on the formation of our Solar System and its uniqueness. 

A fundamental question in this context is the formation of complex abiotic organic 
molecules in protoplanetary disks, such as CH4, formic acid and amino acids. From 
cometary observations in our own Solar System, we know that a small fraction of the 
icy material consists in fact of such organic molecules. During a period of enhanced 
dynamic activity, the "Late Heavy Bombardment", evaporating and impacting comets 
may have delivered large quantities of water and organic material to Earth. But the young 
protoplanetary disks may have been a rather hostile environment for complex molecules. 
The protoplanetary disk is impinged by strong UV and X-ray radiation (stellar activity), 
and a stellar wind with highly energetic particles is driven into the disk at later stages. In 
addition, a massive protoplanetary disk, marginally unstable, will develop spiral waves 
that can lead to local shocks, in which density and temperature can be many times higher 
than in the unperturbed gas. Such shock waves have for example been suggested as a 
viable mechanism to form chondrules Desch and Connolly [1]. In the following, we 
will study the impact of a recurring shock from such a spiral wave on the gas phase 
chemistry. 
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FIGURE 1. Density (a) and temperature (b) structure across the shock. We use the jump conditions 
outlined in Eq.([T]i and (|2]i to compute the peak density and temperature of the gas and apply a fixed 
cooling timescale of 10 minutes. 

THE MODEL 

In a typical region of planet formation, a shock arises from spiral waves that travel at 
a fixed pattern speed as opposed to the speed of the gas i^. As a result there can be a 
large difference in orbital velocity between this spiral pattern and gas moving on nearly 
Keplerian orbits. 

We pick two points at a distance of 5 AU inside a typical T Tauri disk model dsj] 
and fix the initial temperature and density of the unperturbed gas to r(l) = 2094 K, 
n(l) = 10^ cm-3 , Tuv(l) = 2.55 and T{2) = 128 K, n{2) = 10^ cm^^ , tuv(2) = 26.5. 
Within hours after the shock, temperature and density fall back to their equilibrium 
values. 



Chemistry 

The time-dependent chemistry is solved using the variable-coefficient ordinary differ- 
ential equation solver "dvode" [4]. The chemical network is composed of reactions from 
the UMIST database [5] that involve 48 species [3], among which are water and CH4. 
The chemistry has been tested against equilibrium results from Kamp and DuUemond 
risj] for aprotoplanetary disk as well as for a test case of a dark cloud from Millar and 
Herbst M- 



Shock 

We model a local gas volume element at a distance of 5 AU from the central solar-type 
star. The change in gas properties across the shock is defined by a set of jump conditions 



(7+1)^2 



(1) 
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FIGURE 2. Time-dependent abundance of CH4 for the second point in the disk model: the left panel 
shows the evolution until chemical equilibrium is reached, the right panel shows the change in abundance 
during the shock event. 
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where pi and Ti are the unperturbed gas density and temperature, and p2 and T2 the 
physical gas conditions at the shock front. M is the Mach number and 7 the ratio 
of specific heats of the gas before and behind the shock. The latter depends on the 
chemical composition of the gas itself. Fig. [T] shows an example of a shock with a 
cooling timescale of 10 minutes. At a distance of 5 AU the relative velocity between 
the gas and the shock is ~ 10 km/s (taken from simulations of Boss and Durisen L2J). 



TIMESCALES 

Our time dependent chemistry model shows that the chemistry settles back to equilib- 
rium on timescales smaller than the shock orbital period. Fig. [21 and [3] show the time- 
evolution of methane and water abundance at high optical depth (second point in the 
disk model) towards equilibrium and during the first passage through the shock. 



FORMATION OF MOLECULES IN SHOCKS 

We discuss here only the two molecules methane and water. Both of them are formed 
in the shock, reaching abundances up to one thousand (methane) and a million (water) 
times higher than their equilibrium values under the unperturbed gas conditions. 

Analysis of the models reveals that methane is predominantly formed through col- 
lisions between CH3 and molecular hydrogen or water. The main destruction channels 
are collisions with atomic hydrogen, but also to a minor extent photodissociation into 
CH2 and H2. Water is formed via collisions between OH and molecular hydrogen, but 
also through electron recombination of H3O+. Destruction is either through collisions 
with atomic hydrogen (so the direct backreaction that returns the OH to the gas phase) 
or through ion-molecule chemistry, namely H2O -1- Si+ — > SiOH+ -1- H. 
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FIGURE 3. Time-dependent abundance of H2O for the second point in the disk model: the left panel 
shows the evolution until chemical equilibrium is reached, the right panel shows the change in abundance 
during the shock event. 



The efficient formation of water during the shocks has an important impact on the 
cooling timescale. This has already been noted by Desch et al. The presence of large 
amounts of water might lead to very fast cooling timescales of the gas behind the shock 
(only minutes) and this would make the condition for e.g. chondrule formation during 
such shock events very unfavourable. A next step would be of course the self-consistent 
determination of the gas temperature from a detailed heating/cooling balance of the gas. 
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